Обобщение байесовской оценки при учете неоднородности объединяемых данных

При оценке исследуемых показателей байесовским методом предпо­лагалась статистическая однородность априорной и эксперименталь­ной информации. Однако отличие условий испытаний на различных этапах экспериментальной отработки изделий может привести к неод­нородности собираемой статистической информации. Байесовский метод может быть легко обобщен на этот случай [89]. Для такого обоб­щения рассмотрим две гипотезы:

• однородности Uс вероятностью P(U);

• неоднородности U с вероятностью Р(й) = 1 — Р (и).

Тогда в соответствии с формулой полной вероятности апостериорная безусловная плотность вероятности запишется в виде:

P(R/X) = P(R/X, U)P(U) + P(R/X, U)P(U).

Если выполняется гипотеза однородности U, то плотность веро­ятности P(R/x, V) есть байесовская апостериорная плотность вероят­ности. Если выполняется гипотеза неоднородностей U, то целесооб­разно вообще отказаться от использования априорной информации и оценку показателя производить только по экспериментальным дан­ным.

Из выражения для безусловной апостериорной плотности легко получить формулы для оценки и ее дисперсии:

R = R6P(U) + R3KC[-P(U)],

где Дкс — оценка максимального правдоподобия; Д — байесовская оценка;

ДЛ] = Д Д ]P2{U)+Д Л, кс ][1 — P(U)]2.

Из приведенного соотношения можно также найти выражение для доверительных границ, например для нижней доверительной гра­ницы:

J P(R/x) dR = P{U)P(R/x, U)dR + [1 — P(U)]f P(R/ x, U)dR 0 0 0 Проиллюстрируем метод учета неоднородности объединяемых данных на конкретном примере оценки вероятности успешного вы— полнения системой поставленной перед ней задачи. В этом случае, используя условие приемки системы RH = R3 и переходя от интегра­лов к суммам, согласно правилу аппроксимации ^-распределения биномиальным распределением получим выражение для доверитель­ной вероятности у, с которой подтверждается гипотеза R > R^:

image397

где £ (" + «О)- Л, У =1 — уо и

г=0 ‘ ‘

1(;)^_га-^)г = 1-Уэкс

г=0′ ‘

— соответственно выражения для доверительных вероятностей, с ко­торыми подтверждается заданное значение показателя только по экс­периментальным данным и по объединенным данным.

При отсутствии отказов в предварительных и последующих испы­таниях (d = d0 = 0) это выражение приобретает вид:

Р(U)J$+n» + [1 — P(U)]R£ = 1 — у

или Л" = (1 — у)/(1 — P(U) + P(U)( 1 — Yo)) = 1 ■- Уэкв,

где Yo —■ доверительная вероятность, с которой подтверждается за­данное значение показателя А, при проведении предварительных испытаний; уэкв — доверительная вероятность, с которой достаточно подтвердить заданное значение показателя R^ при проведении после­дующих испытаний для того, чтобы обеспечить подтверждение Rj с доверительной вероятностью у по объединенным данным предвари­тельных и последующих испытаний.

Поскольку знаменатель приведенного выражения всегда меньше единицы, то значение уЭкв всегДа меньше у. Следовательно, учет априорной информации всегда позволяет сократить число заключи­тельных испытаний. В табл. 14.3 приведены результаты расчетов не­обходимого числа заключительных испытаний, необходимых для под-

Таблица 14.3

Число испытании, необходимых для принятия гипотезы R > Л,

при различных значениях P(U), у0

т

Число испытаний при = 0,9 и у0, равном

0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

0,95

0,99

0,1

229

228

227

226

225

224

223

222

221

220

219

218

0,2

229

227

225

223

221

218

217

215

211

209

208

208

0,3

229

226

223

220

217

213

209

205

202

198

196

194

0,4

229

225

221

217

211

207

202

1%

190

185

181

179

0,5

229

224

218

213

207

201

193

186

178

170

165

161

0,6

229

223

217

205

196

186

175

184

152

152

145

139

0,7

229

222

213

202

193

185

174

162

147

130

120

112

0,8

229

221

211

200

191

178

164

147

128

103

87

73

0,9

229

219

209

198

190

177

160

145

109

64

37

9

0,95

229

219

208

194

185

169

152

87

37

0

0

0,99

229

218

207

194

180

165

145

112

73

8

0

0

0,999

229

218

207

193

179

161

139

109

69

1

0

0

т

Число испытаний при =

0,99 и у0, равном

0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

0,95

0,99

0,1

298

297

2%

295

294

293

292

291

290

289

288

287

0,2

298

296

294

292

290

287

286

284

280

278

277

277

0,3

298

295

292

289

286

282

278

274

271

267

265

263

0,4

298

294

290

286

280

276

271

263

259

254

250

248

0,5

298

293

287

282

276

270

262

255

247

239

234

230

0,6

298

292

286

274

265

266

254

233

221

221

214

208

0,7

298

291

282

271

262

254

244

231

216

100

189

181

0,8

298

290

280

271

259

247

233

216

197

172

156

142

0,9

298

288

278

267

254

238

221

197

192

133

106

78

0,95

298

288

277

263

249

234

214

181

172

106

68

9

0,99

298

287

276

263

248

230

208

179

156

77

18

0

0,999

298

287

276

262

247

229

207

178

142

70

2

0

тверждения заданного значения показателя при различной степе­ни неоднородности априорных и экспериментальных данных.

При практическом использовании данного метода вероятность Р( U) неизвестна и оценивается как нормированный критический уро­вень значимости, на котором принимается гипотеза однородности объединяемых данных, т. е. P(U) = (а^ — ан)/(ав — ан). Нижним значением уровня значимости ан является рекомендуемый при про­верке статистических гипотез уровень 0,05-0,1. Верхний уровень для двусторонних решающих правил составляет 0,5; для односторонних ав =1.

В заключение следует отметить, что рассмотренный метод при определении вероятности P(U) на основе обработки статистических априорных и экспериментальных данных не исключает полностью проблему смещения объединенных оценок.